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Abstract 

The paper continues a series of publications devoted to the 3D nonhnear locahzed 
coherent structures on the surface of vertically falling liquid films. The work is primar- 
ily focussed on experimental investigations. We study: (i) instabilities and transitions 
leading to 3D coherent structures; (ii) characteristics of these structures. Some non- 
stationary effects are also studied numerically. Our experimental results, as well as the 
results of other investigators, are in a good agreement with our theoretical and numerical 
predictions. 
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1 Introduction 

A Motivation 

Localized nonlinear coherent structures play an important role in many phenomena con- 
nected with the chaotic motion of liquid, in particular, in some regimes of falling films. As 
we discussed in Part I the analysis of experiments reveals the existence of several wave 
regimes, depending on the inlet Reynolds number (Re), (0, (i?e^^^)), {{Re^^'^), {Re^'^'^)) and 
((i?e(2)), (i?e(3))). For water, {Re^^^) ^3-5, {Re^^^) ^ 40 - 70 and (Re^^^) ^ 400. In 
the first interval the instability is too weak for any manifestation to be seen and the 
surface is flat, in the second interval the waves are two-dimensional (2D) and in the 
third interval the surface is covered with three-dimensional (3D) localized coherent struc- 
tures, A-solitons. These 3D solitons are stable and robust and are involved in continuous 
chaotic movement, interacting with each other as "quasi-particles" . This is the so-called 
regime of "surface turbulence" (ST). Note that describing ST in terms of the motion of 
"quasi-particles" results in a dramatic reduction in the number of degrees of freedom. 

In our theoretical investigation we systematically apply the simplified system 
of equations with respect to local film thickness and flow rates in the direction of the 
flow and in the transverse direction, the Kapitsa-Shkadov system. For large Kapitsa 
number 7 = ap^^u^^^^g^^^^ (where a is surface tension, p and u are the liquid density 
and viscosity and g is the acceleration of gravity) the flow is described by the universal 
parameter S (the modified Reynolds number) and the results can be re-calculated for any 
liquid using the relation Re = 7.5l7'^/^^5^/^^. In particular, for water. Re = 66. 15^^^^. 
As Re — )■ 0, the Kapitsa-Shkadov system collapses into the Kuramoto-Sivashinsky (KS) 
equation. 

In Part I - we have examined theoretically the instability of 2D waves with respect to 
a disturbances in the transverse direction. We have also established the 2D-3D transition 
scenarios resulting in fully developed 3D waves whose complex interaction leads eventually 
to ST. In Part II ^, we considered 3D coherent structures of the ST as stationary running 
3D solitary pulses, referred to as A-solitons. We have found these nonlinear solutions and 
have studied their linear stability. It has been found that for the universal parameter 
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in the range S < 0.054 the A-sohtons are unstable with respect to the continuous part 
of spectrum ("radiation instability") and, at 5 > 0.51, A-solitons do not exist. Hence 
their region of existence/stability is 0.054 < 6 < 0.51 which for water translates to 
6 < i?e < 38. 

Among numerous works devoted to the experimental investigation of wave dynamics 
in falling films, only works by Park and Nosoko — and Alekseenko et al. ^ ^ are devoted 
to transitions to the regime of A-solitons and to the properties of these structures. Note 
that for films on an inclined plane with a small inclination angle, 3D solitary waves have 
not been observed (see papers by Liu et al.^ ). 

In the work El water was the working liquid. The authors used artificial perturbations 
in the transverse direction to destroy 2D solitary waves, the critical Reynolds number for 
the transition to 3D waves was evaluated as (Re) ~ 40. At this Reynolds number, the 
most dangerous length of transverse perturbations was about 2 cm. In this work, it was 
stated that disintegration of 2D solitons leads to 3D solitary waves. The work contains 
many nice photographs of the wavy surface, but, unfortunately, quantitative comparison 
with theory is rather difficult: i) a strong influence from the wave interaction complicates 
the analysis of the data, ii) the authors were using the inlet Reynolds number (-Re), while 
for the comparison we need the substrate Reynolds number Re. However, there is a 
qualitative correspondence between the experimental data and our theory. 

As a matter of fact, there are serious difficulties with trying to single out an individual 
"particle" of ST and its parameters from experiments. The main such difficulties are: 
a) the "particles" are the result of a long and complex evolution downstream; b) in the 
ST regime, the "particles" are continuously interacting with each other and it is very 
difficult if not impossible to measure the parameters of an individual particle; sometimes 
it is difficult even to identify them in a randomly disturbed surface. 

In the experiments of Alekseenko et al. SI El, the evolution of 3D localized signals 
was considered for substrate Reynolds numbers Re = 1.25 — 4.7 in an alcohol- water 
mixture with 7 = 404 so that the interval of the universal modified Reynolds number is 
6 = 0.015 — 0.076. In these studies, the authors excited the 3D waves by a short-duration 
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impact with a thin jet of working fluid in the upper part of the film fiow. The major 
trends of 3D locahzed signals were presented and the existence of stationary running 3D 
solitary waves was proved. For the time being, only these experiments contain enough 
information to be compared with the theoretical prediction. 

The main aim of the present work is an experimental investigation of some 3D phe- 
nomena in vertically falling films and a validation of the model developed in Parts I and II. 
The range of Reynolds numbers (-Re) we wish to investigate is 5 to 100 (smaller Reynolds 
numbers were investigated by Alekseenko et al.S^^). First, we study the 2D-3D tran- 
sitions using method ^ and perturbing 2D solitary waves in the transverse direction. 
Second, we artificially create 3D solitary waves by depositing a single drop on the undis- 
turbed liquid surface. A critical aspect of the experiment is that, in contrast to the works 
SI 0, the mass of the deposited drop corresponds to the mass of 3D solitary wave obtained 
theoretically in Part II. To complete the picture of evolution, some nonstationary effects 
are studied numerically. The results of our experiments along with the experiments of 
Alekseenko et al. SI g^^-g compared with our theory in Parts I and II. 

B Experimental methodology 

We are going to investigate the 2D-3D transition and A-solitons for inlet Reynolds num- 
bers {Re) from 5 to 100. For this purpose, we choose a planar vertical channel which is 
25 cm long and 15 cm wide and made of a special mirrored glass. To measure the wave 
characteristics, we choose the fluorescent imaging method by Liu and GoUub ^. The 
method was also successfully used by Vlachogiannis and Bontozoglou 

In order to investigate the transition from the regime of 2D solitary waves to A- 
solitons (2D-3D transition), we applied the ideas of Stini and ^: regular 2D solitary 
waves are created by low-frequency flow rate pulsations; these waves are perturbed in the 
transverse direction by a system of special needles separated by a certain distance. The 
spatio-temporal evolution of such waves are considered in order to flnd conditions of 3D 
instability and transition to 3D solitary waves. 

We devise a novel direct method that allows us to create individual A-structures in 
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our experiments. The basic idea is to deposit drops having the sohton's mass onto the 
hquid substrate. The mass of the sohtary wave is taken from our theoretical calculations 
^. After the drop touches the surface, it is subsequently involved in film movement and 
is rapidly transformed into a stationary running solitary wave. 

C Outline 

In Sec. 2, we obtain the connection between the inlet or pump Reynolds number, (-Re), 
frequently used in experiments, and the Reynolds number based on the soliton's substrate. 
Re, used in our theoretical investigation The ratio (Re) /Re depends on the properties 
of the (3D or 2D) solitary wave, namely its volume and velocity, as well as on the "density" 
of solitons or, equivalently, the number of solitons per unit area. In the regime of ST, 
this ratio is about 10. 

In Sec. 3, we describe our simple experimental setup and a method to measure wave 
characteristics. The planar vertical channel is made of a special mirrored glass and 
distilled water is used as working liquid. Experimental methods to investigate transitions 
leading to 3D solitary waves and to create regular 3D waves on the film surface are 
discussed. The numerical method used to study the evolution of localized signals is 
formulated. 

In Sec. 4, we present our experimental results and compare them with the theoretical 
ones. First, some visual observations of natural waves and of their transitions are pre- 
sented. Second, normal 2D solitary waves, created by periodic pulsations of flow rate, are 
considered. They are excited by special needles in the transverse direction. The evolution 
downstream of these perturbed 2D solitary waves can result in their disintegration and 
further formation of 3D solitary waves. The most dangerous wave length is found. Third, 
in order to create an individual 3D soliton, a drop of liquid with a mass corresponding to 
the equilibrium soliton is deposited directly onto the liquid substrate. Downstream, the 
drop spreads out and takes the form of a stationary running A-soliton. We also perform 
computations to model the creation of the A solitary wave from a localized initial signal, 
or "drop". Good agreement between the theory in Parts I and II and our experiments 
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along with Alekseenko et al. SI S is demonstrated. Finally, we consider a numerical 
solution of simple wave interactions. 

Sec. 5 contains the main conclusions of the work and a discussion of unsolved problems 
for future investigation. 

2 Connection between inlet and substrate Reynolds 
numbers 

The natural control parameter in experiments is the pump discharge flow rate. The 
associated dimensionless parameter is the inlet Reynolds number, (-Re), which is the 
pump discharge per unit channel width divided by the kinematic viscosity of the liquid. 
This Reynolds number is used in the majority of experimental studies of falling films. 
For fixed liquids, it is the only free parameter. 

On the other hand, when we consider individual localized structures, it is convenient 
to introduce another Reynolds number, the Reynolds number of the soliton substrate. 
Re: as we first pointed out in Sec. II of Part I and Sec. Ill of Part I, far from its hump a 
solitary wave decays rapidly to a fiat flow, the so-called "substrate" or "sublayer" . Re is 
based on the substrate's thickness. Note that the substrate Reynolds number has been 
used in the past by several authors, see for example 

Here, we are interested in the connection between the two Reynolds numbers, (-Re) 
and -Re. In the case of wave-less flow, the two numbers are the same, {Re)=Re, but for 
a wavy regime the situation is quite different. In this case, the connection between the 
two Reynolds numbers depends on the number of 2D or 3D solitary waves per unit area 
in the (x, z)-plane or, equivalently, the "density" of these "quasi-particles" . 

Let us first define the volumes of 2D and 3D solitary waves, respectively, 

/oo 
ih-l)dx (1) 
'OO 

and 

/•oo roo 

J^D= I I (h- l)dxdz, (2) 



oo J — oo 



where (2) was first introduced in Eq. (19) of Part II but is rewritten here for clarity. As 
we pointed out in Sec. II D of Part II, even though at large 6 the amplitude of 3D solitons 
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tends to a constant, their volume J^d continues to increase. For 2D solitons, we have 
similar behaviour, i.e. for large 5 their amplitude saturates but their volume continues 
to increase with 5 due to the widening of their humps. In fact, for large 5, the volume 
J2D of 2D solitons is a linear function of 5, much like the volume of 3D solitons, and can 
be well approximated by the relationship 

J2D = 2085 + 0.535 

For completeness, we also give the volume of 3D solitons for large 5 from Sec. II D in 
Part II: 

J^D = 2034(^ - 81.54 

Let us assume that we have a wave-less substrate of unitary thickness with Reynolds 
number Re (the corresponding dimensional thickness is — see Sec. II of Part I). Let 
us also assume that a group of 3D localized coherent structures, A-waves, with average 
distances between them and L^, in the streamwise and spanwise directions respectively, 
slides with a speed c along this substrate. We also define the wave numbers, a — 27: / 
and 13 — 2t:/Lz in the x- and 2;-directions respectively. Let us single out one "quasi- 
particle" and integrate the mass balance equation (5c) in Part I in the frame moving 
with speed c along a period in the 2;-direction and use the fact that the z-component of 
the flow rate, p, is zero far from the hump, 

d c r+^^z^ 

{Q-T- / ^d^) = (3) 



dx J-L,/2 

where Q — Q{x) is a averaged with respect to the z flow rate. 
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Q = — / qdz 



Integrating (3) once with respect to x gives 

g - — / {h- l)dz = 1 (4) 

J-iz J-L^/2 

where the integration constant was determined from the condition that the soliton decays 
rapidly far from its hump to plane-parallel flow h — l,q — 1 (Eq. (7), Part I). Further, 
we integrate (4) once with respect to x: 



1 /•+L^/2 



(g) = 1 + -^-r I I (h- l)dxdz (5) 

where 

{Q) = ^ 

is an averaged flow rate along the x-direction. If and Lz are sufficiently large, the 
double integral on the righthand side of (5) can be approximated by J^d in (2), giving 
the final relationship 

The first term is the contribution from the slowly moving thin substrate; the second 
term is the contribution from the rapidly moving large hump. It depends on the "density" 
of pulses afi, velocity of the wave c and its volume J^d- Mass conservation between the 
flat inlet region and the region of well-developed waves further downstream implies that 
(Q) is equal to the dimensionless inlet flow rate. Hence, the inlet Reynolds number 
{Re), based on inlet flow rate, and the sublayer Re based on sublayer flow rate (whose 
dimensionless value = 1), are connected by the relation 

{Re) = {Q)Re. (7) 

We can also derive for the 2D wave regime an expression similar to (5): 

(g) = i + ^ = i + ^ (8) 

From Eqs. (5) and (6), in the 3D-soliton regime, if either volume J^d ^ or "density" 
a/3 ^ (equivalently "volume" J2D — > or "density" a — > in the 2D-soliton regime), 
then {Q) — > 1 and {Re) Re, as would be the case for a single soliton in an infinite 
domain. However, usually the contribution of the second tern in (5) is much larger than 
the first term and {Re) is much larger than Re. 

The boundary of the 2D-3D transition in experiments is about {Re) fa 40 — 70 for 
water (see Fig. 1, Part I). We should compare these values with our theoretical critical 
(5* — 0.048 or i?e* = 5.5 for water (Sec. IV, Part I), which correspond to instability of 2D 
solitons. In order to do this, we should recalculate Re as {Re) and apply relations (5)-(7). 
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However, this can only be done qualitatively because the "density" of the soliton gas and 
hence (Q) is unknown. Nevertheless, we can use the results of previous experiments to 
estimate the second term in (5) and hence evaluate (Q). According to Chu and Dukler 
the contribution of the second term can be ten to twenty times larger than that of 
the first one, so the majority of the inlet flow rate is carried by solitons. For our range 
of Reynolds numbers, we take the second term to be ten times the first one (see Fig. 12 
of Chu and Dukler's paper). This gives us the rough estimate 

(Re) = lORe. (9) 

Hence, our transition Reynolds number Re = 5.5 is equivalent to (-Re) = 55, which is in 
good agreement with the experimental value {Re^"^^) ~ 40-70 where the 2D-3D transition 
takes place (Fig.l, Part I). 

An important related question is which type of solitons, 3D or 2D, carries more liquid. 
This question, in particular, has important consequences for heat/mass transfer. As we 
demonstrated in Figs. 11 and 12 of Part II, both the amplitude and speed of a 3D soliton 
are smaller than those of a 2D soliton of the same 6. On the other hand, as 6 increases, 
a 3D soliton becomes bulkier, and 3D solitons can potentially carry more liquid than 2D 
solitions for the same 6. By equating (8) to (6), we obtain the limiting value 

/3. = 27r^^^. (10) 

0!3D CsD JsD 

If /3 > 3D solitons carry more liquid than 2D ones; if /3 < 2D solitons carry more 
liquid than 3D ones. Analysis of experimental data shows that the average separation 
distance for 3D waves is much smaller than for 2D, in other words, 3D waves are 
"better packed" in the same domain than 2D waves. As a result a2D / «3d is small which 
together with the fact that J2D/J3D is small (Eq. (1) implies from (10) that P^, is small 
as well). Hence, provided that the average separation distance Lz of 3D waves is not very 
large, and the experiments show that it is not, 3D waves carry more liquid than 2D ones. 
Hence, both cases are possible, depending on experimental conditions. 
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3 Experiment 

Experimental setup. 

We are going to to investigate the range of Reynolds numbers (Re) from 5 to 100. For 
such Reynolds numbers, previous works have analyzed the flow of liquid films falling down 
vertical tubes with diameters from 2 to 4 cm and with water as test liquid. However, 
the typical soliton size in our setup is 3 to 7 cm in the spanwise direction, so such tube 
diameters are too small. We employ an experimental apparatus with a planar vertical 
channel, 25 cm long and 15 cm wide and made of a special mirrored glass, see Fig. [l] 

A digital pump is located in the lower part of the setup since the flow is sensitive 
to external vibrations. Liquid is pumped through a rubber pipe to a supply tank of 
size 25 cm x 15 cm x 10 cm. In the upper part of the tank, there is a feeding device of 
overflow type and the liquid flows out of it through a special polymer net onto the vertical 
test surface. The net serves to decrease the tangential velocity of the liquid. From the 
test section, the liquid flows to the discharge tank. The setup forms a closed loop, such 
that the liquid returns from the discharge tank to the pump. 

Distilled water is used as working fluid. The temperature is kept at 15 ± 0.5°C, so 
that the viscosity is z/ = 1.14 x 10^^ m^/s, the ratio of surface tension to density is 
a/ p = 74 X 10~^ m^/s^ and the Kapitsa number is 7 = 2900. The Reynolds number at 
the inlet, (Re), is measured as the discharge flow per unit width divided by the viscosity 
u. The discharge volumetric flow rate of the pump is varied from 7 to 150 cm'^/s and 
hence the corresponding range of Reynolds numbers is 4.5 < (-Re) < 100. 

It is well known that the main difficulty with an overflow feed in falling film experi- 
ments is the formation of a uniform fiow at the inlet and further downstream. With water 
as working liquid, the film can be very thin. As a result, a small disturbance of the over- 
flow at the inlet leads to the liquid breaking up into drops and rivulets. To remedy this 
effect, we take the following measures: a) the upper edge of the channel is aligned strictly 
perpendicular to the gravity field; b) the channel has a small width, 15 cm, as mentioned 
earlier; c) a polymer net is installed at the inlet to decrease the tangential component 
of velocity; d) in the inlet region, we engrave ten small parallel ditches streamwise, each 
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15 mm long and of triangular cross-section. These measures significantly minimize the 
formation of drops and rivulets. Still, at small Reynolds numbers, (-Re) < 3 — 4, we are 
not able to create a plane-parallel substrate. Finally, polished steel side walls are used to 
minimize the side-wall effect associated with flows in rectangular channels. 

Methods to measure wave characteristics. 
To measure 3D wave profiles and their characteristics, we use the fluorescent imaging 
method by Liu and Gollub ■ The method was also successfully used by Vlachogiannis 
and Bontozoglou El. A fluorescent dye with a low concentration (about 100 - 150 ppm), 
with chemical formula C2oHiQNa20^, is added to the fluid. Such a small concentration 
does not affect the liquid physical properties Eight ultraviolet lamps are arranged 
over the channel surface to illuminate it from above and a chamber is arranged on the 
opposite side of the channel. The absolute accuracy of the method is 8 to 10 /im, hence, 
the largest possible relative error varies from 6 to 8 %. 

The local film thickness can be obtained from the signal by careful calibration. To 
perform the calibration of film thickness versus light intensity, a part of the layer near 
the distributor is used. It remains flat up to the inception point which varied from 7 
to 25 cm. In the measurements, we avoid a small portion of the channel just after the 
inlet to let the flow establish a steady semi-parabolic profile (around ten to twenty film 
thicknesses). The flat portion is used for calibration along with the known discharge of 
the pump. The obtained dependence is shown in Table 1. Between these points, a linear 
interpolation is used, 

/(x, z,t) = K * /o(x, z)h{x, z, t) + Ii{x, z) 

where K is constant, and /o(x, z) and Ii{x^z) depend on the location because of the 
non-uniformity of the ultraviolet light field. 

The digital analogue of the measured physical signal of the surface is recorded. Pho- 
tographs are taken with a high-resolution Panasonic NV-GS75 digital camera. The data 
are recorded on a digital carrier and further processed using a computer with a spatial 
resolution of 720 x 540. In order to measure the wave speed, a series of frames of the 
wave are captured at 1/30 — t/i of a second apart. 
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Typical examples of shadow images of waves, their fluorescence images and cross- 
sections for forced 2D periodic waves with (-Re) = 30 and frequency / = 15.0 Hz and for 
a 3D signal with Re = 6.8 are shown in Fig. 4. 

To measure the wave's amplitude by fluorescence imaging, we use a simple contact 
method in which a steel needle is positioned in the lower part of the channel. The needle 
could move back and forth by means of an adjusting micrometre screw. The flat surface 
of the substrate is used as a reference point. Both electrical and visual methods are 
employed. In the first case, the contact is registered by closing the electrical circuit 
between the needle and contacting liquid. A very small amount of sodium chloride which 
does not affect the hydrodynamical properties of the working liquid is added to it to make 
it conductive. The circuit includes a 4-volt battery and an ampere-meter. To measure 
the maximum amplitude of the wave, we need only one contact with the surface after 
which the experiment is stopped. In the visual method, the circuit is switched off and on 
by the top of the needle where a small coloured particle is positioned. When the surface 
is touched by the needle, the particle is taken by the flow. A microscope is employed to 
attach the particle and to decide if the particle is taken by the flow or not. 

Methods to investigate 2D- 3D transition and 3D solitary waves. 
In order to form 2D solitary waves, localized periodic pulsations in flow rate are super- 
imposed on the main flow 1211121. xhe resultant "artiflcial" waves suppress the natural 
perturbations and a rather regular train of 2D solitary waves travelling downstream is 
generated. In order to induce a 3D instability and the process of disintegration of 2D 
solitary waves into 3D localized coherent structures, we impose artiflcial perturbations 

onto the 2D solitons. For this purpose special needles are placed in the transverse direc- 

~ 3 
tion at the inlet, with a separation distance / from each other, as was done in ^. We 

perform experiments with / from 10 mm to 35 mm in steps of 2.5 mm. 

Let us now discuss how to create 3D solitary waves. Previous falling fllm experiments 

show that 3D solitons can be seen after a long and complex downstream evolution 

(see also Fig. 1, Part I). The key idea of our work is to bypass all intermediate stages 

of evolution and deposit a single drop with a flxed mass onto the smooth fllm surface. 
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The advantage of our method is that we deposit a drop with a mass close to that of the 
theoretical 3D solitary wave. We refer to this method as the "raindrop" method. It was 
prompted literally by the rain on a car windscreen during a traffic jam. This method 
provides a simple, easy to implement and efficient way to successfully create individual 
3D solitons and to measure their parameters as a function of the substrate Reynolds 
number. 

The depositing of drops is done with a special small Finn pipette with adjustable 
volume. We have found this device to be the most appropriate for our purposes. To 
cover the range of soliton volumes obtained in Part II, reproduced here for the ease of 
presentation in Table 2, we employ three Finn pipettes with the following drop volumes: 
0.5 - 10 yuL in increments of 0.1 /iL; 5-40 fiL in increments of 0.5 fiL and 40 - 200 fiL 
in increments of 1 fiL. The smallest diameter of the liquid drop is 2 mm and the largest 
is 7 mm. The calculated dependence of soliton mass on Reynolds number in Table 2 is 
used to estimate the mass of the liquid drop. More specifically, the mass of the liquid 
drop is taken as 1.1 to 1.2 times that of the equilibrium mass. Having a liquid drop mass 
close to the steady 3D soliton mass results in rapid formation of a 3D soliton. 

In the experiment, the pipette is attached to the upper part of the channel such that 
its distance and angle to the surface can be changed. The most painstaking procedure is 
to adjust the pipette in such a way that the moment the required drop size is achieved 
the drop touches the film surface and is subsequently involved in film movement. A 
schematic picture is given in Fig. [7j We have two degrees of freedom to approach the film 
surface, distance and angle to the surface, which are adjusted independently for all the 
fiow rates used. During adjustment, the image of the drop is magnified by a strong lens 
located near the pipette. The distance and angle to the surface are changed by employing 
adjusting micrometre screws. Finally, a table of angles and distances as a function of fiow 
rate is obtained. After the drop touches the surface, it is subsequently involved in film 
movement and is rapidly transformed into a stationary localized wave. 

The important condition in this method is to have a fiat film surface for the deposited 
drop. Recall that the film remains fiat from the inlet up to the inception point. The 
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length of this region is varied from 7 to 25 cm. The drop is deposited onto this part of 
the flow. 

We found that it takes from 5 to 10 cm for the A-wave to reach a steady state — 
the unstable ripples at {Re) < 15 are too weak to affect the formation of the A-wave in 
the upper part of the channel. In this region of Reynolds numbers, the upper part of 
the flow can be assumed to be practically flat. At {Re) ~ 20, room disturbances have 
a sufficiently large growth rate to turn into natural waves; although they continuously 
disturb the configuration of a A-wave sliding on a flat substrate, they are not able to 
destroy it: the A-structure also grows with increasing Reynolds number and is sufficiently 
robust to withstand the natural waves. 

The details regarding the apparatus and the experimental techniques can be found 
elsewhere 

Numerical experiment. 
To complete our investigation, we perform numerical modelling of the non-stationary 
process of A-soliton creation from a drop as done in the physical experiments. As in 
1^^, we use the following Kapitsa-Shkadov system of three nonlinear PDEs: 



dt 5 dx h 5dz h 56 dx ' 



dt^5dxh^5dzh 56^ dz ^ ^ 



dh dq dp ^ 

dt dx dz ' 

where = d"^ jdx^ + d"^ jdx^ is the Laplace operator on the (x, z) plane. All the values 
are dimensionless and incorporate the liquid density p, thickness and averaged velocity 
Mo = (s'/^'o) / (3z^) on the flat substrate. The further transformation, 

^1/3 



32/9i?e2/9 

allows to absorb two parameters, the substrate Reynolds number Re and Kapitsa number 
7 = (Tp~^z/~^/^5f^/^, into a single parameter. 



37/95^1/3- 
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Here, a is the surface tension, p and u are the hquid density and viscosity, and g is the 
acceleration of gravity. 

The numerical method of ^ is used, which is an extension to 3D of the 2D scheme 
developed in the studies of . For a localized 3D signal, the calculations are performed 
in a frame moving with the velocity of the signal. The computational domain is 2Lx x 
2Lz. A sufficiently large size for the domain allows the final A- wave solution to be free 
from "side-wall" effects, or, equivalently, the localized wave has sufficient space to decay 
almost to zero at the boundaries of the domain. We impose boundary conditions at both 
ends of the domain, 

X = —Lx : h = q = 1, 
X = — \- : h = q = 1, p = 0, 
z = ±Lz : h{x,-Lz,t) = h{x,Lz,t), (12) 
q{x, -Lz,t) = q{x,Lz,t), 
p{x, -Lz,t) = p{x,Lz,t), 
The initial conditions are taken in the form, 

t = 0: h = l + Ae'^^''''+''\ q = l, p = (13) 

which simulates a drop on the film surface with parameters A and b describing the 
amplitude and spreading of the drop, respectively. 

4 Results 

Visual observations of natural waves and their transitions. 

First, we consider experiments with waves originated from natural room disturbances. 
A visual examination of the surface shows that it is absolutely fiat at (Re) < 4 - 5, 
confirming previous observations, (Re^^^) 4 - 5. Note that even at (-Re) = 7, the 
manifestation of waves is very weak, see Fig. |3]^a). For larger Reynolds numbers, the 
surface is covered with waves. At any (-Re), experiments show the presence of a fiat 
region near the inlet up to the inception point. At (-Re) > 7, the points of wave inception 
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are located on a line at the following distance from the inlet: (-Re) = 7-25 cm; 10 - 
15 cm; 15 - 10 cm; 20-7 cm. This distance is a little larger than cited in other works 
(see ^) and ^). We attribute the difference to the small width of the channel and 
wall effects. 

For Reynonds numbers {Re) from 7 to 15 - 20, the waves are two-dimensional and 
periodic. For Reynonds numbers from 15 - 20 up to 40 - 50, the waves maintain the basic 
structure of 2D solitons but are disturbed with some 3D modulation, either stationary 
or pulsating. Localized 3D waves can be seen in this range of Reynolds numbers but 
they have small amplitudes and, what is more important, they have a short lifetime and 
disappear rapidly downstream. 

The 2D-3D transition occurs approximately at {Re^"^^) ^ 40 - 50, see Fig. |3](b). The 
3D modulation now is not stationary any more, but grows downstream with eventual 
destruction of the 2D wave and appearance of steady 3D solitons. We note that for 
{Re) > {Re^"^^), the coalescence of neighbouring 2D waves increases the 3D modulation 
and accelerates the destruction process. A typical well-developed 3D regime is shown 
in Fig. ^c). The following interesting visual observation is made regarding the ST: if 
we look at a fixed point on the channel, thus taking the Eulerian viewpoint, the evolu- 
tion appears to be purely chaotic. But if we follow the waves with an appropriate speed, 
thus taking the Lagrangian viewpoint, we can decipher interacting deterministic localized 
structures, which look like "quasi-particles" . In other words, the ST can be considered to 
be a chaotic interaction of quasi-particles which are strongly internally coupled and con- 
tinuously interact with their neighbours. An example of the chaotic wave motion in this 
region is shown in Fig.[6]for {Re) = 60. Several typical patterns can be identified: 1 — 2D 
soliton with 3D modulation prior to being destroyed downstream; 2 — Well-developed 3D 
solitons interacting through their overlapping tails; 3 — Complex coalescence of several 
3D solitons. 

The 3D localized coherent structures slide "rapidly" on a "slowly" moving thin sub- 
layer with Reynolds number Re. In the 3D wave regime, most of the liquid is concentrated 
in the 3D structures and only a small fraction of it is carried by the thin substrate. As 
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we pointed out in the previous section, for our range of Reynolds numbers we take this 
fraction as ~ 10 . As a consequence, the above range of inlet Reynolds numbers for 3D 
wave regimes translates into 4 < Re < 40. 

Forced waves and comparison with our theory. 
The rather irregular picture of natural waves can be improved by imposing artificial 
periodic pulsations in flow rate on the main flow . Further downstream, the pulsations 
are convected into regular 2D periodic waves, at low pulsation frequency the periodic 
waves turn into 2D solitary waves. The initially 2D waves are disturbed in the transverse 
direction by a system of needles placed with a separation / and the regular structure of 
waves is immediately violated. At some distance from the inlet, 3D modulations appear to 
be superimposed on our 2D solitons. These modulations grow downstream and eventually 
destroy the 2D solitons, and 3D localized structures appear at some distance L from the 
inlet. A typical snapshot of such an evolution is shown in Fig. |4| (Re) = 10, I = 25 mm. 

We perform experiments with different /, from 10 mm to 35 mm, in steps of 2.5 mm. 
For sufficiently short separation between the needles, more specifically, I < 10 mm, there 
is no sustained reaction and the 3D perturbations decay downstream. At larger /, wave 
disintegration into localized 3D coherent structures at a distance L is observed. From 
one experimental run to another, this distance varies with an amplitude of approximately 
0.5 cm. Our data are averaged over 20 to 40 points to obtain L. The results of our 
observations for L are presented as a function of the distance / for (-Re) = 10 in Fig.jsj The 
dependence has a pronounced minimum at / = 2 cm which is in good agreement with our 
theoretical prediction ^ of the most dangerous wavelength, 1^ = 1.42i?e^/^ = 1.83 cm, 
and with experimental observations - Moreover, the downstream distance for complete 
disintegration, L, is about 10 to 15 cm, which also fits the theoretical prediction 8 — 15 
cm (Part I Fig. 10). 

In the raindrop method, a 3D coherent structure is formed by depositing a liquid drop 
onto the film surface. The drop is subsequently involved in film movement. We emphasize 
that the drop's evolution results in a stationary 3D solitary wave only if the drop's mass 
is close enough to the equilibrium one. In this case, the drop is rapidly transformed 
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into a stationary localized wave. It takes about 5 to 10 cm from the location where 
the drop touches the film surface for a stationary structure to form. For a larger mass, 
approximately as much as twice the equilibrium mass, a two-hump 3D solitary wave is 
formed. This wave is unstable and decomposes further downstream into two one-hump 
3D solitary waves. If the mass of the deposited drop is much smaller than the equilibrium 
mass, that is approximately one half, then it never evolves into a stationary wave but 
expands in a 3D wave packet. 

A comparison of experimental and theoretical cross sections of A-solitons is given 
in Fig. |8j There is a rather good correspondence between the wave profiles, except 
for the capillary ripple region in front of the wave — theoretical oscillations are more 
pronounced. This can be explained by insufficient accuracy in our measurement in the 
region of capillary ripples, 8-10 fim. Also, the theoretical trough behind the hump is 
deeper than the experimental one, the difference becoming larger at large S. 

The A-wave profiles and other wave characteristics can also be found in experimental 
works HIIS]. In parameters for two solitary waves are presented for 7 = 404: Re = 
2.5, 6 = 0.035, and Re = 3.9, 6 = 0.061. According to our theory in Part II, at 
6 < 0.054 the A-soliton is unstable against the radiation mode, in other words at small 
6, localized disturbances of the plane-parallel part of the A-wave destroy the hump. 
Indeed, the first experimental wave looks unsteady and shows some indication of further 
disintegration. In Fig. |9| we compare our theoretical (a) and experimental (b) wave 
profiles and their cross-sections for the second case. Re = 3.9, 6 = 0.061, 7 = 404. The 
experimental and theoretical velocities are 209 and 213 mm/s, respectively. For this 6, 
the wave is stable and we can find good quantitative agreement between our theory and 
experiment. Our theory also predicts the existence of two-hump 3D solitons and 
their characteristics. Such solitons were found experimentally in The comparison of 



theoretical and experimental data in Fig. [TO] gives reasonably good correspondence. 

Experimental data for the wave speed is given in Table 3. For the first two Reynolds 
numbers, we are not able to create a 3D soliton because of its instability. Instead of a 
nonlinearly coupled localized structure, we observe a spot of expanding chaotic waves 
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confirming that aX S < 0.054 {Re < 6) 3D solitons are unstable, in agreement with 
our theoretical predictions. We are not able to sustain a stationary 3D solitary wave 
at 6 < 0.054 because of a strong primary instability of the fiat region. However, at 
larger 6, the radiation is convected away from the 3D soliton and it is stable and robust. 



The dimensional values from Table 3 are plotted in Fig. [TT] along with the theoretical 
ones. The agreement is really good especially taking into account the simple way our 
experiment is performed. 



Fig. 12 compares the experimentally observed wave amplitudes as a function of 
Reynolds number with the theoretical ones. Again, the agreement is rather good. 

Numerical experiment — comparison with physical results. 
Typical numerical results of the time-dependent evolution for 6 = 0.0908 (for water 



Re = 9.28) are shown in Fig. 13 A localized 'drop' with the initial form described 
by (13) and with a mass 20% larger than the equilibrium 3D soliton mass is placed 
on the flat film substrate at t = 0. At t = 0.05 the drop has a tendency to split 
into two humps; however, provided this event is avoided, instead of disintegration one 
nonlinear hump of the soliton is formed. (The tendency to split into two humps cannot be 
arrested if we take a drop with an initial mass at least twice the equilibrium 3D soliton 
mass.) At time t = 1 the signal acquires capillary ripples (this moment is skipped in 



Fig. 13), and at t = 5 the trough behind the head develops. At t = 12 — 20 the 3D 
soliton is practically formed. One can clearly see at this stage the excess mass draining 
slowly to the back, leaving the equilibrium 3D soliton behind it. The final structure is 
characterized by capillary oscillations/ripples at its front along with moustache-like or 
leg-like structures at the back and a long hollow trough between the legs. The ditch is 
caused by a Bernoulli depression in the moustaches and in fact the numerical experiments 
performed here indicate that it develops simultaneously with the development of the 
moustaches. A completely stationary A-wave running at speed c ~ 4.6 is formed at 
t = 20. 



The results of these calculations are compared with our experiments in Fig. 14 , where 
the wave amplitude hmax is presented as a function of the distance downstream. The 
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theoretical time is recalculated as the distance downstream using the signal velocity 
and physical properties of water. The dimensionless time of 12 to 20 corresponds to 
a dimensional distance of 6 to 10 cm which is in a good agreement with the physical 
experiments. 

The calculations also confirm other experimental results: If the mass of the "drop" is 
at least twice as large as the equilibrium 3D soliton's mass, the drop disintegrates into 
two or several solitons or into a two-hump 3D soliton. If the drop's mass is smaller than 
the equilibrium 3D soliton's mass, the structure sucks up liquid from the substrate as it 
evolves and it takes two to three times longer to reach the final A-wave's steady state. 
Moreover, starting with a small initial mass, the evolution never results in a stationary 
3D wave but in an expanding and growing wave packet. 

The evolution of different shapes of initial signal was also investigated. The constants 
A and b in (13) are varied in such a way that the volume of the initial signal remains 
the same. In all cases, the main stages of the evolution appear to be qualitatively the 
same. Interestingly, the time to reach a 3D stationary structure appears to be almost 
independent of the form of the initial signal. Furthermore, changing the initial mass of 
the drop to 1.5 times that of the equilibrium 3D solitary wave also gives qualitatively the 
same stages of evolution and finally results in a stationary 3D pulse with nearly the same 
time of saturation t ^ 10 to 20. 

It is now well established that the interaction of 2D solitons with smaller objects in 
front of them consists of a coalescence process in which the solitons absorb the mass of 
the smaller objects. This results in an increase in the speed and amplitude of the solitons 
by an amount proportional to the mass of the smaller objects. This coalescence process 
also takes place between 2D solitons of different amplitudes and speeds: the faster soliton 
catches up with the smaller one and absorbs its mass. The coalescence event between 
2D solitons has been observed experimentally (^) and has been described theoretically 

m 



Fig. 15 presents some computational results of such an interaction, but for 3D solitons. 
At it = we place in front of a stationary 3D soliton a localized signal with a mass of 
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I the mass of the sohton. At t = 5, the sohton absorbs the mass of the locahzed signal 
and retains it until t = 25, i.e. for a time interval of At = 20. Recalculation of this 
dimensionless time into a dimensional distance gives roughly 10 cm. The coalescence 
event causes an instantaneous acceleration and increase in speed/amplitude proportional 
to the absorbed mass. The 3D soliton is now an "excited" soliton, i.e. a soliton of larger 
amplitude than that corresponding to the given S. For t > 25, drainage of the excited 
soliton begins and the excess mass leaves the wave. It is expelled from the structure 
through its "moustaches" and eventually takes the form of two small 3D solitons {t = 41). 
These solitons suck up liquid from the substrate and eventually grow into equilibrium 3D 
solitons. This is the key mechanism for generating the staggered "checkerboard" pattern 
seen in many 3D falling film experiments. The beginning of this process is seen at t = 45 
and 49. 

5 Conclusions 

Through simple experiments, we have examined the instabilities and transitions leading 
to the regime of ST. We have shown that if 2D solitary pulses are unstable to transverse 
perturbations, their disintegration results in the formation of 3D A-solitons. For water, 
at sufficiently large Reynolds numbers the most dangerous length is 2 cm; this result is 
close to that obtained theoretically ffl. 

We have developed a new experimental technique, the "rain-drop" method, to inves- 
tigate A-solitons and their characteristics. Wave shapes, amplitudes and velocities have 
been obtained for a wide range of Reynolds number. Numerical experiments modelling 
"rain-drop" formation have also been performed. We have established good agreement 
between the theoretical and experimental characteristics and parameters of A-solitons. 

An important next step in our investigation is the development of a statistical theory 
of the ST. We also plan an experiment and a direct numerical simulation of the 3D wave 
regime in a vertically falling film to model in detail the downstream evolution of natural 
waves, their interaction and statistics. In order to extract the deterministic features 
of the 3D wave regime, the methods of "spatial averaging" and "proper orthogonal" 



21 



decomposition used successfully in ordinary turbulence ^ might also prove very useful 
in the falling film problem. We hope to examine these and related problems in future 
studies. 
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Table 1 Calibration of thickness vs. intensity for the flat film. 
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Table 2 Stretching coefficients k, Reynolds number of substrate and volume of 
equilibrium 3D- "drops", dimensionless and dimensional, as a function of S for water, 
7 = 2900. 
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Table 3 Pipette drop volume J, Reynolds number Re, wave velocity c, and maximum 
thickness of wave hmax, as a function of supply tank volumetric discharge rate, Q (water). 
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Figure 1: Experimental setup. 
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Figure 2: Examples of shadow images of waves, fluorescence images and cross sections: 
(a)Forced 2D periodic waves for Re = 30 and frequency / = 15.0 Hz. (b) 3D-soliton for 
Re = 7. 
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Figure 3: Shadow images of waves generated by room disturbances. The bar is 2 cm long, 
(a) (-Re) = 7, the flow is practically flat, (b)(i?e) = 41, 2D-3D transition and (c)(i?e) = 98. 



30 



6 4 2 

z(cnn) 



Figure 4: Shadow image of periodically excited waves for Re = 10. The waves are 
disturbed in the spanwise direction with needles at the inlet with separation / = 25 mm. 
At point 1, the wave is practically two-dimensional; at point 2 at a distance of 6 cm from 
the inlet, one can clearly see the 3D modulation, at point 3 at a distance L = 11 cm from 
the inlet, the original waves are completely destroyed. 
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Figure 5: The locus of the 2D-3D transition for L as a function of the length of the 
transverse perturbation / for Re = 10. There is a pronounced minimum at / = 2 cm. 
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Figure 6: Typical example of the chaotic wave motion observed for {Re) = 60. This 
section is located between 14 cm and 24 cm from the inlet. There are several typical 
wave patterns: 1 - 2D soliton with 3D modulation ready to be destroyed downstream; 
2 - Well-developed 3D solitons interacting through their tails; 3 - Complex collision of 
several 3D solitons. 
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Figure 7: "Rain-drop" method of exciting 3D soliton with the mass of hquid in the pipette 
corresponding to the equihbrium 3D sohton mass. 1 - Test glass channel; 2 - Pipette with 
adjustable drop-volume; 3 - Successive instants of drop evolution, namely, 3a - drop is 
gradually increasing in size and finally touches the liquid surface and becomes involved in 
film flow, 3b - drop spreads over and turns into an evolving localized signal downstream 
which eventually 3c - becomes a stationary running A wave; 4 - liquid layer flowing down 
the channel. 
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Figure 8: Experimental (1) and theoretical (2) cross sections, (a)i?e = 7, (b)i?e = 8.8, 
{c)Re = 9.7, ld)Re = 13.2, (e) Re = 15, (f) Re = 17.7. 



35 




Figure 9: Comparison of our theoretical (a) and experimental (b) (see Alekseenko et al. 
[5], Fig. 4) profiles of A solitons and their cross sections for Re = 3.9, 6 = 0.061, 7 = 404. 
Experimental wave velocity is 209 mm/s and theoretical velocity is 213 mm/ s. 
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Figure 10: Comparison of our theoretical (a) and experimental (b)(see Alekseenko et 
al. 0], Fig. 3) profiles of two-hump 3D solitons for Re = 2.2, 6 = 0.03, 7 = 404. 
Experimental wave velocity is 102 mm/s and theoretical velocity is 113 mm/s. 
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Figure 11: Phase velocities of A-solitons as a function of Re. The diamonds represent 
our experimental results. 
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Figure 12: Maximum amplitudes of A-solitons as a function of Re. The diamonds 
represent our experimental results. 
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Figure 13: Formation of a A soliton from a localized signal, or "drop", for 6 = 0.0908 
at different times. At t = 0.05 the drop has a tendency to disintegrate. At t = 1.0 the 
capillary ripples form (not shown). At t = 5.0 the trough behind the head forms. At 
t = 12.0 we have reached a nearly steady 3D coherent structure. Finally, at t = 18.0 the 
3D structure takes the shape of a A-soliton. 
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Figure 14: Evolution of the dimensionless amplitude h^ax downsteam, 6 = 0.0908, Re = 
9.28 (water). The line stands for the calculations, the diamonds for our experimental 
results. 
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Figure 15: Coalescence of a A- wave with a localized signal of smaller amplitude and mass 
for 6 = 0.15, (a) t = - 3D-soliton with a a localized signal in front of it; (b) t = 25 - the 
mass is absorbed by the 3D-soliton which becomes an accelerated (excited) 3D-soliton; 
it is ready to swallow up an equilibrium soliton; (c) t = 33 - the extra-mass drains to the 
back of the soliton along its "moustaches"; (d) t = 41 - two small A- waves created from 
the extra mass; (e) t = 45 - "checkerboard" pattern starts to form. 
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